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Abstract 

t— 1 ■ 

In the framework of the lubrication approximation, we derive a set of equations 
describing the steady bottom profile of Leidenfrost drops coupled with the vapor 
pressure. This allows to derive scaling laws for the geometry of the concave bub- 
ble encapsulated between the drop and the hot plate under it. The results agree 
with experimental observations in the case of droplets with radii smaller than the 
£NJ . capillary length R c as well as in the case of puddles with radii larger than R c . 

To cite this article: Y. Pomeau, M. Le Berre, F. Celestini, T. Frisch, C. R. 
Mecanique (2012). 

H . Resume 

6 . 

L'effet Leidenfrost : Description des profiles de gouttes et galettes 

Dans le cadre de 1' approximation de lubrification l'effet Leidenfrost est decrit par 
deux equations couplees pour le profil d'une goutte de Leidenfrost et la pression 
de la vapeur. Ce modele permet de trouver les echelles de longueur caracterisant 
la bulle de vapeur encapsulee entre la goutte et la plaque chauffante. Les profils 
numeriques sont en bon accord avec les observations experimentales, tant pour les 
gouttes de rayon inferieur a la longueur capillaire que pour les galettes de grande 
etendue horizontale. 

Pour citer cet article :Y. Pomeau, M. Le Berre, F. Celestini, T. Frisch, C. R. 
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1 Introduction 



It is a great pleasure to write this piece of science for Paul Clavin. Over the 
years he inspired us in many ways. The topic we have chosen mixes ideas 
of non- equilibrium science, of fluid mechanics and thermodynamics. We also 
predict various properties which agree with the experimental observations. We 
hope that Paul will feel that this piece of science is pragmatic enough ! 

The Leidenfrost effect is named after J.G. Leidenfrost (1715-1794) who wrote 
an article [2], in latin, on his observation that liquid droplets do not touch very 
hot surfaces and so survive much longer than normally expected. As explained 
by Tyndall [3] on nineteenth century, the vapor released in the gap between 
the hot plate and the droplet lifts it and cuts direct physical contact with the 
hot plate. This increases the lifetime of the evaporating droplet because of the 
poor heat conductivity of vapor compared to the one of the hot plate. All this 
works if the droplet is not too heavy. We discussed recently [1] scaling laws for 
small Leidenfrost droplets. This showed the remarkable fact that, for radius 

1 /3 

less than (in order of magnitude) Ri = ( ^ ~ X p[ ) , the droplet takes off from 
the hot plate to reach higher and higher elevations as the droplet gets smaller 
and smaller by evaporation, (ST) being the temperature difference between 
the hot plate and the boiling point of the liquid, rj the shear viscosity of the 
vapor, A its heat conductivity, g the acceleration of gravity, L the latent heat, 
and pi,p v the mass density of the liquid and vapor respectively. Typically Ri 
is in the range of a few tens micrometers (19/xm for a drop of water on a hot 

plate at 400°C), that is much smaller than the capillary length R c = , 
a being the surface tension between liquid and vapor. For water R c ~ 2.8mm, 
namely several order of magnitude larger than Ri. Therefore we shall assume 
Ri «C R c - 

The reference [1] dealt with very small droplets of radii of order Ri or smaller, 
typically from 1/im to 30/xm. In this range it was predicted and observed that 
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the height h of the gap below the droplet increases as R decreases, contrary 
to what is usually claimed (h is also the thickness of the film of vapor between 
the droplet and the hot plate). Starting from droplets with radius larger than 
Ri, it was shown that, as they evaporate and when the radius becomes of order 
Ri, then h becomes of order of the horizontal gap extent / and the droplets 
spontaneously take-off from the substrate, becoming too light to stand the 
upward force generated by the pressure due to evaporation. Because of this 
lift-up the lubrication approximation (for the temperature field and the flow 
in the gap between the droplet and the hot plate) breaks down. 

We focus below on larger droplets, big enough to remain close to the hot 
plate, so that the lubrication approximation applies. We derive first a set of 
two coupled equations for the height h(r) of the droplet bottom surface and 
the local pressure p(r). Increasing the values of R, the solutions are found 
for the four following regimes depending on the location of the droplet radius 
with respect to the radii [Ri, Ri, R c , oo], with Ri defined below. Quasi-spherical 
droplets are described in section 3. They exist for a radius larger or of order Ri 
, but much smaller than described in section 4, they correspond to radius 
much larger than Ri , up to about Ri. In this case the pressure in the gap 
is much smaller than Laplace's pressure (we mean by Laplace's pressure the 
pressure drop across the vapor-liquid interface due to surface tension and equal 
to a{l/R\ + I/R2), R1.2 principal radii of curvature of the surface), therefore 
an " uniform approximation" can be used in the gap, and a spherical shape can 
be assumed close to the bottom of the droplet. In section 5 we shall investigate 
larger droplets, of radius larger than Ri, up to about R c . This requires a more 
complex study because an uniform approximation cannot be used in the gap 
which splits into two domains, a trapped bubble and a narrow neck connecting 
the bubble to the outside. Different scaling laws apply in the trapped bubble 
and in the neck, although the lubrication approximation remains correct in 
both domains. Finally we consider in section 7 the case of puddles, with radii 
much larger than R c . Note that the fluid motion and the temperature field in 
the gap are well described by the lubrication approximation for all cases with 
R larger than Ri. 

Compared to recent publications on the same subject [4] and [5], this work 
seems to be the first one giving (original) estimates of the Leidenfrost effect 
as a function of the physical parameters by discussing the joint phenomena of 
evaporation from the droplet and the viscous vapor flow in the gap between 
the droplet and the hot plate. This led us to the introduction of the length 
scale Ri, which is central in our discussion, as well as the other length scale, 
Ri, depending on it and on the capillary length. To the best of our knowledge 
the set of equations (9), (17) and (19), is used for the first time for solving 
this problem. We notice that in another context similar looking droplet shapes 
have recently been described [5]. In the latter case the drops levitate by air 
cushion above a porous mould through which an air stream is forced. They 
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also display a trapped bubble related to the outside by a neck. In our work 
the flow is a consequence of the evaporation and so of temperature gradient in 
the gap, not an imposed quantity as in [5] , leading to equations and solutions 
with scaling different to ours. 



2 Equations 

2. 1 Velocity field, pressure at the interface and temperature field 

The lubrication approximation in the gap relies on the three following ingre- 
dients. 

i) Stokes equations for the flow in the gap. 

This gap extends mostly in the horizontal direction (x, y), and the components 
of the fluid velocity are (u, v,w),w vertical velocity. The boundary conditions 
are u = v = w = for z = 0, the Cartesian equation of the hot plate. The 
other boundary conditions are on the surface of the droplet, at an elevation 
z = h(x,y). In the lubrication limit, this surface is close to horizontal, so that 
the b.c. (boundary conditions) are u — v — for z = h(x,y). The b.c. for 

AT 

w is Stefan condition, written as w = -jjf- where T tZ \ z= h( x ,y) is the derivative 
of the temperature with respect to z on the surface of the droplet, computed 
on the vapor side (hereafter the notation j\ z will be for Stefan condition 
expresses the conservation of energy: the heat flux XT Z normal to the surface 
of the droplet balances the rate of transformation of liquid into vapor times 
the latent heat. 

The Stokes equations read 



r]V 2 u 



P,x 



(1) 



r]V 2 v 



P,y 



(2) 



and 





(3) 



p being the pressure. The velocity field is divergenceless, so that 



u, x + v, y + w jZ = 0. 



(4) 
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ii) Balance of normal forces on the surface. 

There is another equation for the shape of the liquid surface in the gap. Let us 
consider the case of droplets with radius R much smaller than R c . Depending 
if the pressure generated in the gap by the evaporation flow is of order or much 
less than Laplace's pressure ^ in the droplet, the equation for the surface can 
be discarded (section 3) or not (next sections). 

If the gap pressure is much less than 2|, one can assume that the droplet is 
almost spherical. This yields 

h(x,y) = ho + —, (5) 



h is the point on the spherical surface the closest to the hot plate and r = 
\/x 2 + y 2 is the horizontal distance to this point. The parabolic approximation 
for h(r) is valid in the lubrication limit, h <C R. It is derived in the limit 
r <ti R from the Cartesian equation of a circle, (R — h(r) + h ) 2 + r 2 = R 2 . 

If the fluid pressure in the gap is of the same order of magnitude as 
another equation is needed for h(x, y). This equation results from the balance 
of normal forces on the surface of the droplet. Inside the droplet the pressure 
is dominated by Laplace's pressure (recall that we assume that R is much 
smaller than the capillary radius R c ), whereas on the vapor side the normal 
stress is p — r)w jZ . Therefore the balance of normal forces on the surface of 
droplet inside the gap writes 

2a 

— -(p-riw <z )=a(h. xx + h, yy ), (6) 



which becomes in the axis-symmetric case 

JT ~ (P-Vw,z) = ^(rh r ) r . (7) 



Note that this condition is obviously satisfied by the quasi-spherical profile 
(5) if the pressure in the gap, (p — rjw ;Z ), is negligible with respect to 

in) Laplace's equation for the temperature Field. 

It writes V 2 T(x, y, z) — because we neglect the convective part of this flux, 
assuming the Peclet number to be small. This temperature field satisfies two 
boundary conditions: on the hot plate T(z = 0) = T , and T = 7\ on the 
surface of the droplet , namely for z = h(x,y). In the lubrication limit, the 
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solution of Laplace's equation is 

z \ | Ti z 
h(x,y)J 1 h(x,y)' 

Therefore the vertical velocity on the surface of the droplet is w = —k*Y with 
5T = T — Ti (a positive quantity) and k = 



T = T 



2.2 Pressure in the vapor flow 



Using the above relations, let us derive the equation for the pressure of the 
flow in the gap. By integrating the incompressibility condition from z = to 
z = h(x,y) one obtains, 

<< u tX + v >y » +w(z = h) = 0. (8) 



where << u »= J^~ h udz. In Stokes equation, V 2 is dominated by the 
second derivative with respect to z, the shortest length scale in the lubrication 
limit. Therefore, u is close to the Poiseuille value, 

u = —z(z — h), 
2r] V ; ' 

or << u »= —^f^j-- Once put into the equation (8) it gives, 

V 2 • (-V 2 p) + -Lj- = 0, (9) 



where V 2 = e x -^ + e y -^, e x being the unit vector in the x direction. 

The equation (9) is valid for all situations where the lubrication approximation 
applies. Given h(x,y), it can be written as an Euler-Lagrange condition of 
minimization of the (Rayleigh) functional with respect to variations of p, 



V 



Ra 



= J dx J dy 



h 3 



rjkST 



-p 



(10) 



For axis-symmetric geometries equation (9) reads explicitly, 
frh 3 \ V k5T 
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which can be solved by a double integration, 



p(r) = -(12 V k6T) ) / + Po, (12) 

J nh^n) J h[r 2 ) 

where po is an integration constant fixed by the boundary conditions. Note 
that this expression is valid for any h(r) and requires only that the horizontal 
extension of the gap is much larger than its thickness. 

Let us scale out the various physical quantities which have been introduced. 
As seen in the next section a convenient choice is 

' h s = rTr- 1 ' 2 

< r s = {Rhs) 1 ' 2 = i?f /4 i?V4 (13) 



as units for h, r, and p . For quasi-spherical droplets, this choice readily de- 
rives from the balance between the weight of the droplet and the upward 
force generated by the pressure. We shall see later that it is also pertinent 
for the description of disturbed surfaces (droplets with radius smaller than 
R c ), whereas another scaling will be derived for the description of Leidenfrost 
puddles which forms at R > R C) see section 7. 

With such scalings the equation for the flow in the gap read without any 
physical parameter. In the axis-symmetric case it reads, 

0£*), + s = °- (14) 



and equation (12) becomes 



Thanks to this integral solution, one can see that, given h(r), the integration 
constant po is fixed by the condition that p tends to zero as r tends to infinity. 
The constant p is related to h(r) by the expression 

p„=12/4^/^. (16) 
J nh^njJ h{r 2 ) 
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The balance of the vertical forces on the drop leads to an additional relation. 
The weight of the drop Mg has to be compensated by the vertical force gen- 
erated by the evaporative flow, obtained by integration of the pressure over 
the surface of the sphere in the gap, 

oo 

F z = 2n J drrp(r) = Mg. (17) 



Using the relations (13) the balance of vertical forces becomes in dimensionless 
form, 



2 x 



3 



drrp(r). (18) 



o 



2.3 Pressure versus curvature 



From w :Z = one finds that p is larger than the viscous stress r/w jZ , by a 
factor l 2 /h 2 , I horizontal extent of the gap . Therefore, when the lubrication 
approximation applies, the contribution of the viscous stress to the balance of 
vertical forces can be neglected, and equation (6) yields 

\-v = °^\h (19) 



which writes in the axis-symmetric case, 

2(7 1 

— -p = a(h r .2 + -h <r ). (20) 



After it is written with the units given in (13) this equation becomes 

2 - = h r 2 + hl >r = ^(rh r ) tr , (21) 



where all the variables are scaled, making appear the dimensionless number 

e = r 7/2 r; 3/2 r~ 2 , (22) 



which is the ratio of the pressure in the vapor flow to Laplace's pressure in the 
drop. The scaled coupled set of equations (14)-(21) together with appropriate 
boundary conditions make up the model we shall consider henceforth. We 
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recall that they are derived thanks to the scalings (13) appropriate for drops 
radii smaller than R c , as studied in sections 3-4-5, whereas another set of 
scalings will be derived for larger drops (puddles), leading to the same system 
of dimensionless equations, see section 7. 



3 Solution in the case of an undisturbed surface 



In this section we solve the equations for the lubrication limit in the case where 
the droplet remains almost spherical, that implies to neglect the term £p in 
equation (21). This eliminates the need to derive the shape of the surface of 
the droplet in the gap by using the balance of normal forces. The present case 
(undisturbed sphere) is a fairly standard application of lubrication theory. We 
start using the original variables and show that the scalings proposed in (13) 
naturally yield parameterless equations. 

The result of the integration of the right hand side of equation (12) with h(r) 
given by equation (5), is 

12rjkR5T„. r 2 . 

P( r ) =P0 713 G (- 



4/ig v /i iT' 
where G(a) is the numerical function defined as 

" da' ln(l + a'/2) 



f da' in i + c*'/2 



To have a pressure tending to zero at r tending to infinity one must take 
p Q = 12 ^ 5T G(oo), whence the result, 



which yields the vertical force generated by the evaporative flow, F z defined 
in (17), 



3r]kR 2 8T 

By writing that this force balances exactly the weight of the sphere, 

oo 

2tt J drrp(r) = Mg, (24) 
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one finds h = 




that agrees with the scaling proposed in (13) for the 



gap height. Whenever the surface of the droplet is close to a parabolic cylinder, 
equation (5) for h(r) in scaled variables becomes 



Notice that po is now a pure number constrained by the condition that p(r) 
tends to zero as r tends to infinity, and that ho in equation (25) is also a pure 
number defined by the balance of vertical forces on the droplet (18). 

The set of equations solve the quasi-spherical problem if the lubrication ap- 
proximation applies, namely if the height h Q is physically much less than R, 
which requires R ^> R[. It also assumes that the pressure in the gap is neg- 
ligible compared to a/R, because we assumed the relation (5). This requires 
£ < 1, or R < Ri with 



which is about 330/im for water over a plate heated at 400°C . Otherwise, 
one has to determine the shape of the droplet in the gap, namely the function 
h(r), as done next. 



4 Solution in the case of a disturbed surface 

In this case the surface of the droplet in the gap is not a spherical cap because 
the pressure generated by fluid motion there is not negligible with respect to 
Laplace's pressure inside the droplet. Compared to the previous case, we have 
to solve the same equation for p already written in (9), but the profile being 
unknown, we need also to solve equation (6) with the convenient boundary 
conditions. This second equation is derived from the balance of normal forces. 
Assuming the equation (9) solved, the pressure p is known. 

In the range £ ~ 1 the new equation to be considered is the equation for the 
curvature of the droplet in the gap, which should replace the simple relation (5) 
used in section 3 for the range £ 1 characterizing an undisturbed spherical 
droplet. Finally, if £ is not small, equations (9) and (19), or their scaled form 
(14) and (21), make together a pair of equations allowing to obtain the droplet 
profile h(r) and the pressure p(r) in the gap. 

Some properties can be derived without explicitly solving the equations. From 
the integral solution for the pressure given in (15), the pressure is a decreasing 



(25) 



Ri = (RfRl) 1/7 



(26) 



10 



function of r because h(r) is positive. Because it has to tend to zero at r 
infinite, the pressure is positive and decreasing. From equation (21) the mean 
curvature ^(r/i r ) r is an increasing function of r, tending to 2 (in dimensionless 
units) as r tends to infinity. This excludes in particular very large values of 
this curvature at finite values of r in the limit £ large. 

The condition of balance of vertical forces can be transformed into a condition 
for the behavior of h(r) at large r. In the dimensionless version of the equa- 
tions, this condition is derived by integrating both sides of equation (21) from 
zero to a large radius with the element of integration (rdr). The final result 
yields the following condition valid at r tending to infinity: 

yln(r) + .... (27) 



This is a way of expressing that the vertical force on the droplet is equal to 
the uncompensated vertical component of the capillary forces. 

One can reduce the equations to a single one for h(r) with a closed set of b.c. 
By simple algebra, one derives: 



The Laurent expansion of h(r) near r = reads h(r) = h + a\r 2 + a 2 ln(r) + 
a 3 r 3 + ... where h and ax_ 3 are free coefficients, although the coefficients of 
the next order terms in the expansion, like 6r 4 , etc. can be derived order by 
order from equation (28). The coefficients and a 3 must vanish to make the 
solution smooth. The two remaining free parameters ho and ai are fixed by 
the asymptotic behavior of the solution at large r. It reads 

f 

h(r) = cr 2 + d\n(r) H — - + 

where c, d, and / are free parameters. Two parameters are constrained by the 
condition that a 2 and a 3 vanish, but the solution that we shall display has 
0-2 = «3 = and so we shall not consider them anymore, and call a\ simply a. 
The next order terms are also derivable from c and d by order by order solution 
of the equation (28). The b.c. on the shape of the surface imposes c = 1/2 
and the condition for the balance of vertical forces imposes d = -f . This 
yields two conditions for two free parameters, ho and a. The exact solution 
of equations (28) fits well the analytical expansions (27) and (25) (valid for 
r < 1 and r > 1 respectively), as shown in Fig.l where the solid line coincides 
with the dotted line, as distinct from the spherical profile located above. 
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h(r) 
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3 
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Figure 1. (a) Bottom profile of a Leidenfrost drop for £ = 0.18 or R = 200 \im. The 
numerical solution of equation (28) in red solid line melting with the asymptotic 
expansion (27) in black dashed line, are located below the spherical profile (solid 
purple line). Initial conditions h r 2(0) = (2 — £)/2, ho = 1.375. The physical vari- 
ables can be recovered by using the scaling lengths h s = 5.86/um, r s = 34am. (b) 
Decreasing pressure versus axial distance. £p(r) solution of equations (21). 

5 Limit £ large 

The limit £ small is somewhat trivial since in this limit one can neglect the 
term £p in equation (21) which has the simple solution h(r) = h + y, and 
one is back to the case of a spherical droplet of unit radius, as expected. 

An obviously interesting limit is the limit of a large £. In this limit, the solution 
splits into two different domains. Those domains are derived from an analysis 
of the solution of the equations in this limit of a large £. The results presented 
are consistent with the equations, although they cannot be considered as ob- 
vious consequences of them. Between r = and r = r c (to be found) the 
surface of the sphere is like a trapped bubble with a negative curvature and 
where the pressure p is almost constant. The radius r c is such that the trapped 
bubble solution crosses the hot surface, which is obviously impossible. To get 
rid of this crossing, other scalings must be used locally and a neck replaces 
this crossing, as shown below. From equation (21) the mean curvature of the 
surface is also constant, and negative (which is the only possible choice as one 
can check). 



5.1 Unbalanced pressure hypothesis 

Assuming for the moment (something that will be shown not to be correct) 
that this constant value is not 2/£, the curvature of the surface of the trapped 
bubble should of order 1, one finds that h b (r), namely the value of h(r) in the 
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trapped bubble is equal to 



7 / \ 2-(po 2 

hb{r) = r + h - 

4 

We assume now that the main contribution to the upward force is from the 
bubble, something to be checked at the end. This yields from equation (25) 
that 

2= 4h 

This gives a trapped bubble which is an almost perfect spherical cap with 
constant pressure inside. It should be connected to the outside by a neck 
smoothing the solution near r = r c , as necessary because the integral contri- 
bution to the right-hand side of equation (15) scales formally like r 2 //i 4 , of 
order £ -3 in the bubble although p scales like which makes it dominant. 
Therefore inside the trapped bubble the pressure is constant and equal to p . 
The balance of vertical forces yields 

7 Pot -2 

ho = -r. • 

3irp 

This solution does not work however because it cannot be matched with the 
solution in the neck. In other words we have to discard the hypothesis of 
spherical cap (for the surface of the trapped bubble), because the continuity 
of slope of the surface at the transition between the neck and the trapped 
bubble cannot be insured. 



5.2 Description of the trapped vapor bubble 



The only possibility remaining is that the pressure inside the bubble is at 
leading order (with respect to £) equal to 2/£ to balance Laplace's pressure 
inside the droplet, plus a small contribution depending on r and balancing the 
curvature of the droplet there. Therefore we assume that, inside the trapped 
bubble, 

2 

P = 7 + Pb, 

where Pb(r) is much smaller than At leading order p — |. This makes it 
straightforward to derive the radius of the trapped bubble from equation (18) 
which set out the balance of vertical forces. One finds 

rl = f- (29) 



This radius is actually the radial distance between the axis at r = and 
the neck where the pressure makes a transition from its value 2/£ inside the 
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bubble to zero outside of it. Therefore the radius inside the trapped bubble 
is of order ^Z 2 , the order of magnitude of r c . The equation relating p b to the 
curvature of the surface reads 

- iPb = h b , r 2 + -h b>r . (30) 



The other equation relating h b and p b is derived from (14) and reads 

fiR + £=°- (31) 



Defining scaled quantities £~ 1/,2 r, £~ 3 / 5 h b and ^ +1 ^p b makes disappear any 
small or large parameter in the differential equations to be satisfied by those 
quantities. It means also that, inside the trapped bubble, the height is of order 
£ 3 / 5 and that the correction to the leading order constant pressure (2/£ in the 
original variables) is of order £~ 7 / 5 , negligible in the large £ limit, as expected, 
with respect to the leading order contribution 2/£. Note also that even though 
the scaling law h b ~ £ 3 / 5 seems to imply that the lubrication approximation 
does not hold because h seems to be much bigger than r c ~ £ 1,/2 , this is not so 
because physically h and r are originally measured with different unit lengthes. 
We shall come back to this at the end of this section. 

With the scaled quantities (written the same as the original quantities) the 
equations to be satisfied are equation (31) with p b instead of p and 

~Pb = h b:T 2 + -h htr . (32) 



The condition to be satisfied by this set of equations is h b (r c ) = with 

r c = (j^J ^ . Actually this crossing is unphysical. It defines the large distance 
behavior (in inner variables) of a neck solution connecting the trapped bubble 
with the outside. 



5.3 Neck region 



Let 5 = r — r c be the local coordinate in the neck and h n (S) be the local height. 
In the neck the pressure is of order l/£ because it has to tend to p = 2/£ 
on one side (in the trapped bubble) and to zero outside. From equation (18) 
because p is of order h n (r) should scale like 5 2 . Assuming the neck to 
be much less extended than the trapped bubble, namely that 5 r c , one 
finds that, inside the neck h is much smaller than in the bubble, namely much 
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smaller than £ 3 / 5 , and 

h n , fi = 2-&(6). 

This is to be completed by equation (14) relating the pressure and the height 
inside the neck. As far as the order of magnitude with respect to £ is concerned 
this last equation is consistent if p scales like l/£, 6 scales like £ 1//6 and h n (S) 
scales like £ x / 3 (This power law £ x / 3 makes, as expected, the height in the neck 
much smaller than the height of the trapped bubble the latter being of order 
£ 3 / 5 ). Let us introduce local (overlined) quantities by absorbing the scaling 
laws in multiplicative factors, 

p = ip 

< 5 = r 1/6 5 (33) 

h n £ I h n . 

V 



The equations to be satisfied by the overlined quantities are purely numerical 
(namely without large or small parameter) and read 



0, (34) 




and 



2-p 



h T 2. 

n,o 



(35) 



The boundary conditions for the pressure are simple to write. For 5 very large 
negative (that is on the side of the bubble), the scaled pressure should tend 
to po = 2, and as 5 tends to plus infinity (namely outside the bubble) the 
pressure should tend to zero. The asymptotic conditions for h(5) are dealt 
with later. 

Supposing h n known, one can integrate equation (34) to obtain: 

p© = -12/=§-/j| T+ p , (36) 



The integration limit r = of the original problem is pulled to 5 = — oo for 
the overlined variable. The jump of pressure from p = p Q = 2/£ to zero takes 
place almost exclusively across the neck. This yields p = 2 as a constant of 
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integration. The other b.c. for the pressure is p — > as S tends to plus infinity. 
It becomes the following condition for h n : 



+00 — 



2 = p = 12 [4^- [J%. (37) 



The solution in the trapped bubble and in the neck should merge somewhere. 
This merging occurs in a region where the solution of the equations on either 
side have a common power law behavior with respect to the distance to r c . 
This power law behavior is easier to see on the neck side. It has to do with 
the way the solution behaves for 5 large negative. The equations to be solved 
are (34) and (35). Let us introduce p 1 = 2 — p. The equations to be solved 
become 




and 

Pi — h j2. 

They have as an exact solution the power law, 
and 

pi = Mr 6/5 , 

that makes an acceptable solution because, at large negative 5, p l decays to 
zero. Therefore the pressure tends to 2, as it should. 

Outside the neck, the solution should merge with the asymptotic spherical 
droplet in a region where the pressure tends to zero. For the neck solution 
this writes h = 5 2 + a5 k . Introducing this expression in equations (34)-(35), 
we obtain the asymptotic expansion valid for large 5 

K^5 2 + ^-r A + ... (38) 



and 



P 



(39) 
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Figure 2. Experiment and theory for £ = 80. (a) experimental setup ; (b) 2D plot 
of the portion of surface below the drop (interference pattern) for water Leiden- 
frost drops of radius R = 1.14mm on a substrate kept at To = 400° C; (c) ex- 
perimental profile h — ho (circles on the black curve) of the droplet bottom along 
the horizontal axis of curve (b). The red curve inserted below is the theoretical 
solution of axis-symmetrical equations (14)-(21) for £ = 80, with initial conditions 
h r 2(0) = — £~ 2 / 5 , /io = 0.9£ 3 / 5 plotted in physical variables using the scalings 
h s = 2.44/im, r s = 53/im, p s /cr = 0.07/im -1 defined in (13) . 



The above description of a vapor bubble encapsulated below the concave part 
of the droplet in the large £ regime was done under the hypothesis of axis- 
symmetric profile, although experiments report mostly non axis-symmetric 
ones [6]-[7]. Therefore a quantitative comparison between numerical data and 
experiments is out of our scope. Nevertheless we have found a fair agreement 
between numerical solutions of equations (14)-(21) and the experimental pro- 
files for large £, as illustrated in Fig. 2. In this figure the experimental set-up 
is depicted in (a) (see the detail in [6]). A high speed camera is used to image 
the interference fringes between the drop- vapor and substrate- vapor interfaces 
which are visible in (b). The relative height profile below the drop along the 
horizontal of figure (b) is reported in (c), where the experimental asymmetri- 
cal black curve ( with circles) qualitatively agrees with the axis-symmetrical 
solution of equations (14)-(21) obtained for the same value of parameter £. 

As shown in Figures 2(c) the trapped bubble size and height are clearly compa- 
rable in the numerics and in the experiment. Furthermore the trapped bubble 
connects with the outside through a very narrow neck, as predicted above. 
The narrow neck is observed for large £ values corresponding to the domain 
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Figure 3. Numerical solution of equations (14)-(21) for £ = 1000, or R = 2.36mm. 
The height h (red solid curve ) and pressure ap (blue dashed curve) are plotted in 
physical variables, with a = 140 a scale factor. Initial conditions h r 2(0) = —0.15, 

h = e /5 - 

Ri < R < R c as illustrated for £ = 1000 in Fig. 3. In the latter figure we also 
report the pressure profile in the gap (dashed curve) which decreases from 
the center of the bubble to the outside of the neck, where it vanishes. The 
theoretical curves are drawn with the physical variables by scaling back the 
variables h,p help to the expressions written in the next subsection. 



5.4 Large £ domain in physical variables 



The case £ large considered here is such that the droplet radius fulfills the 
condition Ri <C R <C R c - The length scale h s for the height of the droplet 
above the heated plate, and the length scale for the horizontal distances r s are 
given by the relations 13. Therefore the radial extent of the trapped bubble is 

n ~ R\R c y\ (40) 



which is also the radius of the disc of contact at equilibrium of small nonwetting 
droplets. The vertical thickness of the trapped bubble is 

h b ~ # 8 / 5 i?f /5 (i? c )- 6 / 5 , (41) 



The gap elongates in the horizontal direction as R increases, with an aspect 
ratio h b /r b = R~ 2/5 R^ 5 R~ 1/5 becoming: smaller and smaller as R increases. 
Therefore the lubrication approximation remains correct as £ gets bigger and 
bigger, which is equivalent to increase R at constant R c and Ri. Therefore 
the range of applicability of the theory at large £ relying on the lubrication 
approximation extends all over the range [Ri, R c ] and stops to be valid when 
R becomes of the same order of magnitude as the capillary radius. 
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6 Evaporation of droplets in the various regimes 

An obvious application of the ideas presented before is the derivation of the 
lifetime of the droplet. This is related to the evaporation flux from the droplet. 
This flux depends itself on the range of parameters where the radius lies. At 
very small radii one expects that, because the droplet is well above the hot 
surface, the evaporation flux is very close to the one for an isolated droplet 
(without the hot plate) in a vapor at temperature T + ST at large distances. 
In the case of an isolated droplet the total mass evaporation rate is given in 
order of magnitude by J is ~ R 2 p v w with w = -^-T :Z and T z ~ T/R, that 
gives 



This mass loss rate is the one explaining the decay of the square radius pro- 
portional to time, the well-known D 2 law. We shall consider below the mass 
loss of matter out of the droplet in various range of radii. The corresponding 
law of decay of the radius is easy to get in each case and we shall not do the 
calculation. Actually we shall limit ourselves to find whether the main con- 
tribution to the evaporative mass comes from the quasi spherical part of the 
droplet or from the evaporation in the gap. 

Let us compare to Ji s the evaporative mass in the vapor film between the 
droplet and the hot plate in the regime Ri <C R <C Ri. In this regime, the 
area of the film is much less than 47rR 2 , the area of the drop, but the tem- 
perature gradient in the film is much larger than ST/R, so the normal speed 
of evaporation is much bigger than what it is for an isolated droplet with the 
same ST. Therefore it is not obvious which mass loss rate, Ji s or Jfu m is the 
largest. The total mass loss rate in the film is of order l 2 p v w ~ ^tt- From 
the relation I 2 = Rh one derives immediately that Jfu m ~ Jj S where J is is 
given in (42), an unexpected result valid in the regime Ri <C R <C Ri. 

Let us estimate Jfu m in the regime Ri <C R <C R c . We use again the relation 
J film ~ r\p v w ~ r bTK~ together with the estimates r& (radius of the film) 
and hb (thickness of the film or of the trapped bubble), given by equations 
(40)- (41). Inserting those estimated into the expression of Jfu m pertinent for 
this case, we obtain 




(43) 
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/ n \ 7/5 

Therefore, in the range Ri <C R <C R c the ratio Jfum/Jis is of order ( |M , 
much bigger than one, so that the flux is dominated by the contribution from 
the film. We did not consider the flux coming from the neck domain. It is rather 
straightforward to show that it is negligible compared to Jf# m in the limit £ 
large and <C R, equivalent to i?, C C R c - Our theoretical prediction 
is confirmed by the experiment, as shown in Fig.4. In this figure the experi- 
mental data (circles) are obtained from droplets evaporating over a thin brass 
substrate (thin enough to be curved and then stabilize the droplet) healed at 
300°C. The red curve is the best fit to equation (43) in the whole domain 

/ \ 1/5 

Ri < R < R Cl leading to a mass evaporating rate very close to iww) > 

an expression that could be derived from Jfu m ~ 7rr 2 p v w by using the relation 
(29) for r c = r b /r s . 
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Figure 4. The experimental mass evaporation rate (points) are plotted together with 
the theoretical prediction (solid curve). The red curve is the best fit to equation (43), 
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very close to Jfu m 



From this discussion one can conclude that, even though the low heat con- 
ductivity of the vapor lowers the rate of mass loss of the droplet, it does it in 
different ways depending of the range of droplet radius one considers. 

A droplet beginning with a radius of order R c will decay to a (much smaller) 
radius Ri after a time of order t* ~ [RlR]) 1 ' 5 ^^. If one substitutes R 

for (R 3 Rj) 1 / 5 in this estimate one gets back the standard D 2 law. This is 
consistent with the acceleration of the evaporation: to evaporate an isolated 
droplet of radius R c would take a time of order R\ X ^ T , much longer than the 
time t* just estimated (recall that R c ^> Ri). 
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7 Leidenfrost puddles 



At equilibrium, drops of radius larger than R c standing on a flat plate become 
flat puddles. Of course they cannot be characterized anymore by their radius 
only, because their shape is not spherical. The relevant quantity to characterize 
them is their volume. If this volume V were spherical, it would define a radius 

R—\^) ■ Differently from the drop case above, we consider below physical 
variables, in the limit R ^> R c , with R so defined. In this limit, we assume 
that the equilibrium puddle of liquid has a circular shape of radius of order 
r c and thickness of order R c . From the estimate of the volume of the puddle, 

/o3\l/2 

the radius r c = (-^-J can be considered as given and much bigger than R c . 



7. 1 description of the vapor layer 



We consider now the situation of a flat puddle hovering on a hot plate by the 
Leidenfrost effect, and derive the structure and thickness of the vapor layer 
between the plate and the puddle. Like in the case of droplets of radius in the 
range [Ri, R c ] we assume the puddle to be very close to its equilibrium shape 
at the given volume. This has to be checked at the end to yield a solution 
consistent with the assumptions and with the underlying physics. A direct 
consequence of the assumption of closeness to the equilibrium shape is that 
the pressure in the puddle near the bottom is almost constant and just equal to 
Archimedes hydrostatic value pigH, H being the height of the puddle, of order 
R c , so that we shall simply replace H by R c in the coming order of magnitude 
estimates. Therefore the situation is similar to the one studied before in the 
sense that the pressure in the liquid (above the vapor film) is constant, the 
Archimedes's pressure replacing the Laplace's one. Let p ~ PigR c be this 
Archimedes pressure. 

The equation (20) for the pressure reads now 

Po ~P = <r(h r 2 + -h r ), (44) 



As we did before we assume that p is the constant p plus a small (with 
respect to a parameter to be found) part p 1 which balances the curvature 
term in equation (44). Therefore 

-Pi = <T(h r 2 + ^h tr ), (45) 



This is to be completed by the equation (11) for the inhomogeneous part p\{r) 
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of the pressure 



The two equations can be reduced to a dimensionless form by taking 



R c 



h c = (Rfr A c R~ 2 ) 1/5 = (RfR 6 R~ 4 ) 1/5 (47) 
p lc = (R?r^R 7 Y /5 , 



as unit for r, h and for pi/a. The two resulting equations read 

-Pi = (h <r 2 + ^h r ), (48) 



and 



rh 3 
"12 



Pi, 



(49) 



The condition pi <C po is satisfied if r c 3> (RcRi) 1 ^ 2 , which is satisfied because 
we assumed r c >> R c ^> Ri. 

There is another constraint on this solution: the film of vapor must be a 
thin layer underneath the puddle. Therefore its height must be much smaller 
than the height of the puddle. This implies R c 3> h, equivalent to R c 3> 
i?f /5 r^ 5 i?- 2 / 5 , or to R\ > R\r\. This puts an upper bound on the radius of 

/„ \3/4 

the puddle : jf- <C l^J , which is compatible with the conditions R c <C r c 
and Ri <C R c - Lastly a neck makes the transition between the trapped bubble 
and the outside. In the neck the equations to be satisfied are formally the same 
as equations (34) and (35), except that the scaling are derived differently. The 
constant pressure in the trapped bubble is now p Q = pigR c - One obtains 
the same equations as (34) and (35) by taking as units of width, height and 
pressure, 

5=(R c R l ) 1/2 

< h = Ri (50) 

[P 2 
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hfr) Pi(r) 




Figure 5. Solution (a) h(r) and (b) pi/po of equations (48)-(49) for R = 7.5mm. 
Numerical initial conditions h r 2(0) = — 1, ho = 4.1 in scaled variables, (c)-(d) 
Experimental (points) and theoretical (solid lines) values, of the height h in (c) and 
neck radius r ne( ± in (d) versus R, of the vapor bubble below the drop. The piecewise 
theoretical curves are plotted in physical variables, from equations (40)-(41) for 
Ri < R < R c , and relations (47) for R c < R < R ma x, both without any fit. 

The b.c. are the same as before, namely p tends to 2 as 5 tends to minus 
infinity and p tends to as 5 tends to plus infinity. 

It is remarkable that this puddle solution is given by the solution of the same 
set of equations, but with different scaling parameters. Remember that it has 
a limited range of existence because, physically, the wider is the puddle, the 
larger is the pressure in the trapped bubble, so that it reaches the top of the 
puddle at the limit of its domain of existence. It could be that at still larger 
masses of fluid, steady solutions have chimneys of vapor crossing the puddle 
from bottom to top. In reality it is likely that this corresponds to the onset 
of boiling, an unsteady phenomenon in general. Moreover the trapped bubble 
could be Ray leigh- Taylor unstable as soon as r c gets bigger than few R c , an 
instability that could be counteracted by viscosity of the vapor in the gap. 
Therefore steady solutions in this range of values of r c 3> R c are at best only 
indirectly connected to real life phenomena. An example of numerical solution 
for the puddle is given in Figure 5 where curves (a)-(b) display the profile 
and pressure for R = R max , the largest radius of stable water puddles. We 
plot in curves (c)-(d), solid lines, the theoretical values of the height and neck 
radius for the two domains [Ri, R c ] and [R c , R m ax] as derived in paragraphs 
(5) and (7) respectively. The points reported on these figures are experimental 
results corresponding to Fig. 2 of [7], that are in very good agreement with 
our theoretical predictions drawn without any fit to experimental data. Let 
us precise that while the variable r max in [7] equals R for R « R c , it is 
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close to R(j^) 1 / 2 for R ^> R c . In between we have extrapolated the r max (R) 
dependance in order to compare the data of [7] with our theoretical results. 



7.2 evaporation rate 



The mass evaporation rate of the vapor layer below puddles can be derived 
by using similar argument as used for droplets in section 6. It writes Jfu m ~ 
r 2 c p v w ~ r cTif where r c and h c are the the radius and height of the film. 
Introducing the physical values written in equation (47), we obtain 

X5T ( R 9 \ 1/5 (R\ 4/5 (R c \ 3/5 
Jfttm ~ —r- I -5-53 I =J is -5- -5- , (51) 



R c Rf J \RiJ \R% 



which is still much larger than the evaporation rate of the equivalent spherical 
isolated drop. Let us compare this mass loss rate (out of the vapor layer) with 
the rate out of the top of the puddle, and also with the rate out of the neck 
region. To calculate the temperature gradient close to the top surface of the 
puddle, we see the puddle as a disc of radius r c and temperature T—ST inserted 
on a plate heated at temperature T. This leads to a vertical temperature 
gradient of order ST/r c because r c is the unique scale length of this problem. 
It follows that the evaporation rate is J top ~ r cTr~ ~ ~JfUm & t the top 
of the puddle. On the other hand the mass loss rate in the neck region is 
Jneck ~ "^TX ~ ^r^Jfiim- Using equations (47)-(50), these two relations 
lead to the ratios 



Jtop ( Ri 



J film \VRR~c 



3/5 

I , (52) 



and 



Jneck ( RlR c 



J film \ R? 



1/10 

! , (53) 



that are much smaller than unity for R t R c < R < R max ( recall that R ma x 
is of order few R c ). We conclude that the evaporation mostly takes place in 
the vapor layer below the puddle, as it was shown for droplets in the range 
Ri < R < R c in section 6. In summary the evaporation process happens 
through the vapor bubble for any drop or puddle whose radius R belongs to 
the domain [R h R max \. 
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8 Summary and perspective. 



This contribution explains the scaling laws for the Leidenfrost phenomenon, 
sweeping the domain of small droplets to large puddles. Our approach relies on 
a scaling analysis of the fundamental equations in the lubrication approxima- 
tion, which has a wide domain of applicability for explaining the Leidenfrost 
phenomenon. A significant restriction on the validity of our approach is a con- 
straint on the control parameters: we assumed R c ^> Ri, in agreement with the 
data for the experiments done in Nice. It could be however that this inequality 
is not satisfied in other experiments, opening the way to other scaling laws. 
Another very interesting question is the transition to boiling. It has to do with 
the extension of our approach to large puddles. Experimentally this could be 
related to the much studied and still mysterious effect of film boiling. There 
one has to deal with the occurrence of physical contact between the liquid and 
the hot plate by the breaking of the continuity of the vapor film, something 
beyond our approach, as we assumed the hot plate to be at fixed temperature, 
above the boiling temperature of the liquid. Likely the understanding of the 
transition to film boiling requires a solution of the heat transfer equations in 
the solid also. The ratio of heat conductivities of the vapor and the solid could 
be used as a small parameter for this problem. 

Another limit for the applicability to real experiments of the concepts of fluid 
mechanics and heat transfer, as used in the present work, is the Knudsen 
limit: the mean-free path of molecules in ordinary conditions in air is in the 
micrometric range, and could be of the same order or even bigger than some 
length scales of the Leidenfrost phenomenon. Indeed the smallest scale is the 
thickness of the film of vapor underneath the droplet. If this thickness gets 
noticeably smaller than the mean free path, the correct physical picture for 
the flow there is by the direct solution of the equations of kinetic theory, the 
Boltzmann equation in principle or simplified versions of it, like the BGK 
model [9]. In the limit of a mean-free path much bigger than the thickness 
of the vapor layer, one recovers a rather simple description: the molecules 
bounce on both sides of the thin layer and so make a Brownian motion in the 
horizontal direction. This is described in the lubrication approximation by a 
diffusion equation with the time and the horizontal coordinates as variables, 
and the number density per unit horizontal conserved quantity. This 

has to be matched with the regular continuum mechanics picture (Stokes 
equation and Laplace's equation for the temperature field). This transition 
to a rarefied gas situation could explain some of the observation of splashing 
of droplets at low pressure [8]. In this much studied problem of the impact 
of droplets on solid surfaces (in the absence of any temperature effect), the 
pressure in the trapped bubble should depend on the history of the collision, 
and some scaling laws derived in this paper could be valid. For example small 
impacting droplets at low speed could remain quasi spherical and so could be 
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dealt with the same lubrication approximation that we used. At larger impact 
velocities it is likely that a concave trapped bubble could show up too. All 
this will be the subject of future investigations. 
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